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Comparative Physiological and Molecular Analyses of Intraspecific 
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Abstract; Stipa purpurea Griseb, the dominant species of alpine steppe, is widely distributed across the large pre- 
cipitation gradient on the Tibetan Plateau. It is possible that because of local adaptation, S. purpurea populations 
from different habitats along this precipitation gradient respond differently to drought, which may affect their respon- 
ses to climate change. To explore the problem, in the present study, we investigated the physiological and molecular 
response of S. purpurea seedlings from two different populations (Pulan & Cuoqin) to 14-d drought stress and subse- 
quent recovery. The results showed that the relative water content, chlorophyll fluorescence, content of osmoticant 
proline and malondialdehyde (indicator of oxidative stress), accumulation of reactive oxygen species, antioxidant 
enzyme activities and the expression of drought-related genes all changed under drought treatment and went back to 
the control levels in the subsequent recovery in plants from Pulan. However, these patterns were quite different in 
plants from Cuogin, in which these traits changed by inconsistent degrees and did not return to pretreatment levels 
after rewatering. The results demonstrated that the plants from Pulan had greater resistance to drought stress com- 


pared with those from Cuoqin, which had a larger mortality rate ultimately. Combating the differences of offspring in 
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response to drought and the habitat distribution of parents, we considered that genetic basis has been obtained in re- 


sponse to precipitation difference among S. purpurea populations. The results help to understand the adaptation and e- 


volution of S. purpurea to the special environment and the effect of climate change to this botanical system. 


Key words: Stipa purpurea; Tibetan Plateau; Drought; Adaptation and evolution; Climate change 


Stress constantly challenges the survival and de- 
velopment of plants, and plants have in turn evolved 
many adaptations to different environments. Among 
the different stresses, drought stress draws much at- 
tention because it is a major factor that decreases 
plant growth and productivity ( Yordanov et al., 2000; 
Aranjuelo et al., 2011). Plants endure drought stress 
mainly by avoiding tissue dehydration, while main- 
taining as high tissue water potential as possible, or 
by tolerating low tissue water potential. In fact, 
plants initiate an integrative mechanism during this 
process, involving phenotypic plasticity, alterations 
in biochemical enzyme activity or proteomic dynam- 
ics, and stress-responsive changes in gene expres- 
sion ( Nicotra and Davidson, 2010). Under drought 
stress, water loss in plants can be minimized by roll- 
ing leaves (Turgut and Kadioglu, 1998; Terzi et 
al., 2013) , thick leaves (Hu et al., 2012) , devel- 
oped trichomes (Fu et al., 2013), or more rigid 
cell walls ( Moore et al., 2008) etc. In addition, the 
plant can increase the root investment, such as in- 
creasing the root depth ( Moore et al., 2008). Con- 
sistent with phenotypic plasticity, physiological and 
biochemical adjustments are involved in the toler- 
ance to a low tissue water potential (Chaves et al., 
2003). Plants close their stomata to reduce water tran- 
spiration and photosynthesis when exposed to drought 
(Chaves et al., 2003). The osmotic compounds, 
like proline, also regulate osmotic adjustment in plant 
adaptation to drought (Morgan, 1984). A series of 
scavenging mechanisms, including antioxidant en- 
zymes, can restrict the over-accumulation of reactive 
oxygen species (ROS), which can damage cell 
membranes and other cellular components during 
drought stress ( Scandalios, 1993). Correlated with 
both morphological and physiological divergence, 


water shortage is reported to trigger drought-respon- 


sive gene expression, which is characterized by the 
induction of signal transduction ( Yamaguchi-Shi- 
nozaki and Shinozaki, 2004, Huang et al., 2012) , 
channel protein genes (Johanson et al., 2001; Luu 
et al., 2007; Danielson and Johanson, 2008) , tran- 
scription factors (Nelson et al., 2007; Qiu and Yu 
2009; Tripathi et al. 2014), and other defense 
genes (Hare et al., 1999; Catala et al., 2007; Sun 
et al., 2013) , to subsequently maintain homeostasis 
(Fordyce, 2006). However, plant drought stress re- 
sponses are complex biological processes that are 
controlled by multiple trait loci, and most studies to 
date have focused on model and crop plants, neglec- 
ting alpine plants which survive under harsh environ- 
mental stress and are possibly highly valuable for 
their potential resistance genetic resources. Clearly, 
many challenges remain for the comprehensive un- 
derstanding of the mechanism of plant tolerance to 
drought stress. 

Climate warming has caused substantial changes 
in temporal and spatial precipitation patterns in man- 
y regions of the world (Zhuang et al., 2010). The 
Tibetan Plateau, famous as “The third pole” and 
“Roof of the world” , is closely connected with glob- 
al climate change (Carlyle et al., 2014; Wang et 
al., 2014). Although there are plentiful freshwater 
resources in the Tibetan Plateau, and this region 
serves as the headwater for many of Asia’ s rivers 
(Wang et al., 2009), the distribution of the water 
resources is predominantly dependent on its topo- 
graphic configuration and atmospheric circulation 
pattern. The present atmospheric circulation in Asia 
is dominated by a westerly and monsoon wind sys- 
tem, which results in an approximate rainfall gradi- 
ent from east to west on the Tibetan Plateau ( Klein 
et al., 2004; Shen et al., 2008) , with accompanying 


changes in ecosystem types from alpine meadow, 
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steppe to desert (Ni, 2000). Meanwhile, global cli- 
mate change and human activities are potential rea- 
sons for the accelerating drought trend on the Tibet- 
an Plateau. Tree-line data indicates that extreme 
drought conditions have progressed into the north- 
western Tibetan Plateau during the past 300 years 
(Gou et al., 2006, 2007). Therefore, rainfall is 
one of the most important factors determining the 
structure and function of plant populations, commu- 
nities and alpine ecosystems on the Tibetan Plateau, 
indicating that relevant studies should be very help- 
ful for ecosystem management, conservation, and 
development. However, most previous studies have 
focused on macroecological effects ( Bothe et al., 
2010, 2011; Zhu et al., 2011; Yang et al., 2012; 
Sun and Liu, 2013; Wu et al., 2013; Yang et al., 
2013; Wang et al., 2013; Fang et al., 2014), 
whereas very little is known with regard to the effects 
of rainfall on the evolution and adaptation of the par- 
ticular species or populations on the Tibetan Plateau 
at the physiological and molecular levels. 

Stipa purpurea Griseb. is an endemic and domi- 
nant perennial grass species in the alpine steppe and 
meadow environments that are widely distributed in 
the Tibetan Plateau, Pamirs Plateau and high moun- 
tains in Central Asia ( Yue et al., 2011). Because 
of the wide distribution and strong resistance to 
cold, drought and high winds (Zhou et al., 1987) , 
it plays a vital role in safeguarding soil and water re- 
sources, acting as a windbreak and preventing ero- 
sion (Yue et al., 2011). With global climate change, 
S. purpurea is an available candidate to study alpine 
plants response to drought stress and mine the resist- 
ance genetic resources. As S. purpurea has an ex- 
tremely wide range of distributions on the Tibetan 
Plateau, it may be locally adapted to regional envi- 
ronmental variation. In the present study, we investi- 
gated the physiological and molecular divergence of 
plants of S. purpurea from two population response to 
drought stress in the laboratory and analyzed how the 
difference was correlated with their respective pheno- 


type and acclimation to habitat environments. We 


supposed that the S. purpurea population living in the 
more arid regions might have relatively stronger tol- 
erance to the same degree of drought as a result of 
heredity effect. The results should help us to better 
understand how physiological and molecular diver- 
gence contributes to drought adaptation across natu- 


ral populations of S. purpurea on the Tibetan Plateau. 


1 Materials and methods 
1.1 Seed collection and cultivation 

The mature S. purpurea seeds were collected in 
August 2013 at the time of seed release from two 
populations of western Pulan and eastern Cuogqin re- 
spectively on the Tibetan Plateau (Fig. 1). After 
brought back to laboratory, the seeds were dried un- 
der a constant condition of 15 °C , 15% air humidity 
for one week. 

The full seeds of two populations were picked 
out randomly and the seed awns were removed. Then 
the seeds were sown in flowerpots [ diameter = 6. 8- 
9.4 cm (from bottom to top), height =8 cm] with 
equal weight of humus soil, watering to keep the soil 
sufficiently moist. Each pot contained thirty seeds. 
The pots were placed in a greenhouse under 12-h- 
light/12-h-dark , 28/20 °C and a relative humidity of 
40% - 60%. For each population, leaves of ten-pot 
seedlings were collected and grinded for proline, 
MDA, antioxidant enzyme activity measurement and 
RNA extraction. 3 one-pot seedlings were used to 
analyze change of morphology, RWC in leaf and 
chlorophyll fluorescence respectively. Three leaves in 
one pot were selected randomly for H,O, and O,7 
detection. Three independent replications were car- 
ried out for all sample measurements. 

1.2 Drought treatment 

Before drought treatment, the seedlings grew 
for another three weeks watered every day uniformly. 
When the plants reached the 3-leaf seedling stage 
(Liu et al., 2013), the plants were withheld water 
for two weeks until the plants of one population showed 


apparent death. Then the plants were rewatered for 
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Fig. 1 Geographical information and location on the precipitation distribution map of the two seed collection sites 


another two weeks; after recovery, the plant mortali- 
ty for each population was counted. The leaves of 
plants under drought for 0 d and 14 d and rehydra- 
tion for 14 d were sampled for the following measure- 
ment and analysis. 
1.3 Water content measurement 

The relative water content (RWC) of the leaves 
was determined as; RWC =[ Fresh weight (FW) - 
Drought weight (DW) ]/( Turgid weight-DW ). To 
measure turgid weight, leaves were kept in distilled 
water in darkness at 4 °C to minimize respiration los- 
ses until they reached a constant weight ( Rivero et 
al., 2007). DW of the plants was obtained after 48 h 
at 70 C in an air oven (Rivero et al., 2007). 
1.4 Analysis of chlorophyll fluorescence 

Chlorophyll fluorescence was analyzed as previ- 
ously described (Bai et al., 2011), with a pulse- 
amplitude modulation chlorophyll fluorometer ( Heinz 
Walz GmbH, Effeltrich, Germany). Briefly, at the 
time of sampling, S. purpurea seedlings were dark-a- 
dapted for 30 min to measure the maximum quantum 
yield of photosystem II (PSII; F,/F,,), and F/F, 


was determined for each sample by analyzing a whole 


pot plant. The maximum fluorescence (F’,,) was re- 
corded by a 0.8 s pulsed light at 4 000 umol s”' 
m”, and minimal fluorescence ( F,) was recorded 
during the weak measuring pulses. 
1.5 Proline and malondialdehyde content meas- 
urement 

Proline content was measured as previously re- 
ported (Bates et al., 1973), with some modifica- 
tion. Approximately 0. 2 g of fresh leaves was homog- 
enized in 10 mL of 3% aqueous sulphosalicylic acid, 
and the homogenate was centrifuged at 2 000 x g for 
10 min. Then, 2 mL of the extract were reacted with 
2 mL of acidic-ninhydrine and 2 mL of glacial acetic 
acid for 1 h in boiling water; after this, the reaction 
terminated in an ice bath. The reaction mixture was 
extracted with 4 mL toluene, mixed vigorously with a 
test tube stirrer for 15 -20 sec. The chromophore 
containing toluene was aspirated from the aqueous 
phase, warmed to room temperature and the absor- 
bance read at 520 nm using toluene for a blank. The 
proline concentration was determined from a stand- 


ard curve and calculated as mg g' DW. 


The malondialdehyde (MDA) content was de- 
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termined as described previously (Li et al., 2014). 
Approximately 0.5 g of fresh leaves were homoge- 
nized in 10 mL of 10% trichloroacetic acid (TCA) 
and centrifuged at 12 000 x g for 10 min. Then, 2 mL 
of 0. 6% thiobarbituric acid in 10% TCA were added 
to an aliquot of 2 mL of the supernatant. The mixture 
was heated in boiling water for 30 min and then 
quickly cooled in an ice bath. After centrifugation at 
10 000 x g for 10 min, the absorbance of the superna- 
tant at 450, 532, and 600 nm was determined. The 
MDA concentration was expressed as nmol g`' DW. 
1.6 In situ H,O, and O, detection 

The in situ detection of H,O, and O,” were per- 
formed using a previously reported method with some 
modifications ( Able, 2003). The amount of H,O, 
and 0,” was detected with 1 mg mL”! of diaminobenz- 
idine (DAB) and 10° M nitro-blue tetrazolium (NBT) 
respectively. Leaves were vacuum-infiltrated in 10 
mL solution for 2 h and then cleared in boiling etha- 
nol (95%) for 10 min. Then, samples were stored 
and examined in 95% ethanol. 
1.7 Antioxidant enzyme activity assays 

The activities of superoxide dismutase (SOD; 
EC1. 15.1.1), catalase (CAT; EC1. 11.1.6) and 
ascorbate peroxidase ( APX; EC1.11.1.11) were 
determined as previously described ( Nakano and 
Asada, 1981; Jiang and Zhang, 2001). Approxi- 
mately 1 g of leaves collected from treated S. purpur- 
ea seedlings were homogenized in 10 mL extraction 
buffer (50 mmol + L™' sodium phosphate pH 7.0, 1 
mmol: L' EDTA, 1 mmol+ L' DTT, 1 mmol - L” 
GSH, 5 mmol-L™' MgCl, + 6H,0, 1% PVP-40 and 
20% glycerin). The homogenates were centrifuged at 
12 000 x g for 10 min at 4 °C, and total soluble pro- 
tein content of the supernatants was measured by the 
Bradford method ( Barbosa et al., 2009). 
1.8 RNA extraction and RT-PCR analysis 

Total RNA was extracted from the leaves of dif- 
ferent samples using TRIzol reagent ( Invitrogen ) 
and treated with RNase-free DNAse (Takara). The 
RNA concentration was determined using a Nanodrop 


1000 ( Thermo Scientific product, USA). MMLV re- 


verse transcriptase ( Promega ) was used to synthesize 
the cDNAs. The alternatively spliced fragments were 
amplified from the cDNA using primers listed in Ta- 
ble 1, which were designed based on the transcrip- 
tomic sequencing results. The amplification was per- 
formed under conditions as follows; 94 °C , 4 min; 94 
C, 30s; 52-58 C , 30s; 72 °C, 30s; 28-35 cycles; 
72°C, 10 min (Table 1). The amplification products 
were separated by electrophoresis on 2% agarose gels. 
1.9 Data analysis 

Statistical analyses were performed using the 
statistical Software Package for Social Science (SPSS) 
version 18.0. One-way ANOVA for all variables was 
used for testing the treatment differences. Differences 


were considered to be significant at the P < 0.05. 


2 Results 
2.1 The phenotypic difference between plants 
from two sites 

After exposure to drought stress for 14 days under 
the same conditions, the plants of Cuogin showed ob- 
vious wilting compared with those of Pulan (Fig. 2A). 
Some of the plants died after another 14-d rehydra- 
tion with significant divergence between plants of two 
sites. Actually, 16.2% of the seedlings from Pulan 
were dead while mortality rate of Cuoqin reached to 
78.2% (Fig. 2B). The phenotypic results fully dem- 
onstrated that there was difference between the plants 
from two sites in response to drought stress. 
2.2 The changes of physiological status of the 

plants 

The RWC of the plants from both Pulan and Cuo- 
qin had high values (83.9% and 83.3% respectively ) 
before drought treatment and dropped to dissimilarly 
low levels (32.0% and 13.6% respectively ) after 
drought stress for 14 d, with obvious difference be- 
tween two sites (Fig. 2C). After rewatered, the RWC 
of plants from Pulan returned to a relatively high level 
(69. 1%) while that of Cuoqin still had a very low 
value (35.1%) (Fig.2C), indicating that maybe 
the water absorption system of plants from Cuogin 


was destroyed by the stress. 
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Table 1 The information of genes related to drought stress used for RT-PCR 


GenBank 





Gene name A Primer F (5'-3') Primer R (5'-3') T PCR cycles 
accession number /C 

Actin1 KM216249 GCTGGATTCTGGAGATGGTGTC TTACTCATTCACCACTACGGCTG 52 28 and 35 
APX1 KM201379 CACGAGGAAGAATACACCC GCATAAAGCTCCACATAGC 50 30 
APX2 KM201380 CCTTCGGCACCATGAAGTG GCCTCAGCATAGTCAGCAAA 54 30 
CAT KM201378 CACGCCTTCAAGCCCAAC ACGAACCTGTCCTGCCTGTC 56 30 
DREB1 KM201387 GGACGTTCCCAACTGCTC CGGTTCACCTTCTATCGG 52 32 
DREB3 KM201388 GTTGGCTCGTACTAACTTCC ATCATCGGTTCACCTTCTAT 56 32 
PSCS KM201383 TCATAACTGGCGTCATTC TI'GTCCACTCCCTCGTAG 50 30 
PIP1 KM201389 TCCGAGGACAAGGACTACAAGG GGACGAAGGTGCCGATGAT 56 30 
PIP2 KM201390 CATTGACTTCGAGGAGCTGACC GGTGCTTGTACCCGATGACG 56 30 
SOD1 KM201381 CTTATTTGAGCAAGAGGG CTGATGGCATTTCGTGTA 48 30 
SOD2 KM201382 GAAGCACCACGCCACCTA GCTCCCAGACATCAATCCC 54 30 
WRKY4 KM201385 AGCCTTCACTGAGCCTTGACC CCTTCTGCCCGTACTTCCG 58 30 
WRKY11 KM201386 CCCTCCCGTITCTTCTCCTCCT GGGCTTCTGTCCGTACTTGCG 58 30 
WRKY17 KM201384 CAGAAGCACGTTAAGGGA GAGGCTCACGGTTAGGTT 52 30 
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Fig. 2 Changes in plant phenotype and relative water content of leaves during the process of drought treatment and subsequent recovery 


of Stipa purpurea from two sites. A; The change of plant phenotype of Stipa purpurea from two sites; B; The change of plant mortality 


rate of Stipa purpurea from two sites; C; The change in relative water content of leaves in Stipa purpurea from two sites. 


Error bars indicate SE. Means denoted by different letters were significantly different (P< 0.05) (B & C) 


The chlorophyll fluorescence was diversely in- 
fluenced in the course of drought stress and subse- 
quent recovery. The F/F „ of plants from Pulan just 
changed from 0. 86 to 0. 52 during the drought treat- 
ment (Fig.3); nevertheless, the F/F „ reduced from 


0. 86 to 0. 05 in the plants from Cuoqin (Fig. 3). In 
the recovery stage, the F/F „ returned to an average 
value of 0. 81 and 0. 20 respectively in the plants of 
Pulan and Cuoqin (Fig. 3) , suggesting the different 


influence degree under drought stress. 
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Fig. 3 Effects of drought stress and subsequent recovery on leaf photosynthesis in Stipa purpurea from two sites. A; F,/F „ images. 


The pseudocolor code depicted at the bottom of the image ranges from 0 (red) to 1.0 (purple). The experiment was replicated 


three times with similar results. One representative result was shown; B; Average F/F „ values. F,/F „ was determined for 


whole plants. Error bars indicate SE. Means denoted by different letters were significantly different ( P < 0.05) (B) 


2.3 The changes of proline and MDA content 

The proline content of plants from both Pulan 
and Cuoqin kept at a relatively low level before 
drought stress (1.51 and 1.52 mg g' DW respec- 
tively) (Fig. 4A). After treated by drought stress for 
two weeks, to our surprise, the proline content of 
plants from Pulan was drastically induced up to 
10.11 mg g™' while it rose up to 6.07 mg g™ in 
plants from Cuogin, which were 5.7 and 3.0 fold 
respectively more than the corresponding controls 
(Fig. 4A). In the subsequent recovery, the proline 
content in both plants recovered by certain degree, 
but they were still higher than the controls (1.96 
and 2. 68 mg g`' respectively) (Fig. 4A). 

The MDA content was diversely influenced in 
the course of drought stress and subsequent recovery. 
It changed from 4. 18 to 17. 09 nmol g’ DW during 
the drought stress in the plants from Pulan (Fig. 4B) , 
nevertheless, it increased from 4. 07 to 22. 83 nmol 
g DW in the plants from Cuoqin (Fig. 4B). In the 
recovery stage, the MDA content returned to the 
control level in the plants from Pulan (3. 34 nmol g"! 


DW), but it remained at a higher level than the 


control in the plants from Cuogin (12. 86 nmol g`' 
DW) (Fig. 4B). 
2.4 The accumulation of H,O, and O,- 

Both the H,O, and O,` in plants from two sites 
were accumulated under exposure to drought stress 
and dropped to a low level after rewatered, but they 
both showed different change degree between plants 
of two sites (Fig. 4C). Both H,O, and O,” were in- 
duced to a much lower level and returned to be more 
similar with the control in the plants from Pulan 
compared with those from Cuoqin (Fig. 4C). The re- 
sults were consistent with the MDA content, which 
reflected the oxidation degree of cell membrane by 
the excessive ROS accumulation. 

2.5 The dynamic of antioxidant enzyme activities 

The SOD activity, as well as APX activity, 
were induced significantly under drought stress and 
returned to the control level roughly in the plants 
from Pulan (Fig. 5A). However, both SOD and 
APX activities did not increase but decreased slight- 
ly after drought treatment in plants from Cuogin 
(Fig. 5A) ; when the plants werewatered again, the 


activities of both the two antioxidant enzymes contin- 
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ued to reduce to a lower level compared with the 
control (Fig. 5A). The CAT activity also rose greatly 
in plants from Pulan when the plants were treated by 
drought stress and went back to the control level af- 
ter rehydration treatment, while it just increased 
slightly under drought stress and the value reduced 
to the control below after rehydration in the plants 
from Cuogin (Fig. 5A). The results indicated that 
the antioxidant enzyme system in plants from Cuoqin 
was inefficient and might be partly destroyed in the 
course of drought stress. 

2.6 The expression changes of genes related to 

drought stress 


To investigate the expression patterns and func- 





15 
mm Pulan A 

= == Cuogin 
8 d 
‘oo 104 
on 
= 
g c 
B 
= 
5 54 
£ b 
£ a a i 

0 





Control 


Drought Recovery 


H2072 accumulation 


Control Drought 














Recovery 


tions of genes in response to drought, RT-PCR was 
performed for several drought-related genes. The re- 
sults indicated that different genes showed various 
expression changes pattern with two common points 
(Fig. 5B). On one hand, nearly all genes’ expressions 
were induced by different degrees under drought 
treatment in both sites of plants ( Fig. 5B). On the 
other hand, most genes had higher expressions in 
plants from of Pulan compared with those from Cuo- 
qin after exposure to drought stress ( Fig. 5B). More 
specifically, two SOD genes ( SOD1, SOD2) and 
one CAT gene were all induced obviously under 
drought stress and returned to the control levels after 


rewatering in plants from Pulan, but they showed 
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Fig.4 Accumulation of proline, malondialdehyde (MDA) and reactive oxygen species (ROS) ( H,O, and 0,” ) in Stipa purpurea from 


two sites during the process of drought treatment and subsequent recovery. A; Proline content changes in Stipa purpurea from two sites; 


B: MDA content changes in Stipa purpurea from two sites; C; Jn situ detection of changes in leaf H,O, and O,” levels of Stipa purpurea 


from two sites. Error bars indicate SE. Means denoted by different letters were significantly different ( P < 0.05) (A and B) 
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slight induction under drought stress and dropped 
below the control levels after watering again in plants 
from Cuogin (Fig. 5B). Two APX genes ( APX1, 
APX2) were also induced under drought stress and 
went back to the control levels in plants from Pulan; 
on the contrary, expression of the two genes did not 
increase but decrease after drought treatment in 
plants from Cuogin (Fig. 5B). The delta-1-pyrroline- 
5-carboxylate synthase-like ( P5CS) gene showed 
obviously high expression in both plants from Pulan 
and Cuogin (Fig. 5B). Three WRKY family genes 
(WRKY4, WRKYII and WRKYI7) had similar ex- 
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pression changes in the plants from the same site; 
nevertheless, they showed better induction and re- 
covery in plants from Pulan than those of Cuoqin 
(Fig. 5B ). Interestingly, the similar results were 
found in two dehydration responsive element bin- 
ding protein (DREB) family genes ( DREB1 and 
DREB3) (Fig. 5B). Two aquaporin genes ( PIPI 
and PIP2), which belong to the class of plasma 
membrane intrinsic proteins ( PIPs), also showed 
higher expressions in plants from Pulan compared 


with those from Cuoqin after exposure to drought 


stress ( Fig. 5B). 
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Fig.5 Changes in antioxidant enzyme activities and drought-related genes expression in Stipa purpurea from two sites during the 


process of drought treatment and subsequent recovery. A; The changes of antioxidant enzyme activities in Stipa purpurea from 


two sites; B; The changes of drought-related genes expression in Stipa purpurea from two sites. Actin was included as a cDNA 


loading control. The abbreviations Co, Dr and Re at the bottom represent Control, Drought and Recovery respectively. 


Error bars indicate SE. Means denoted by different letters were significantly different ( P < 0. 05) (A) 
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3 Discussion 

Intraspecific variation in drought response has 
been observed from a lot of species under different 
growth conditions ( Walker and Miller, 1986; Sand- 
quist and Ehleringer, 1997, 1998; Beikircher and 
Mayr, 2009; Lu et al., 2009; Hamanishi et al., 
2010). These findings show the effects of environ- 
mental heterogeneity on the modification and accli- 
mation in plant characteristics. As S. purpurea is dis- 
tributed along a representative precipitation gradi- 
ent, the characteristics of response to water shortage 
or excess also may exhibit significant differences in 
different populations, with a hypothesis that plants 
from more arid area may possess greater drought tol- 
erance. In the present study, the results of morpho- 
logical and physiological status in the course of 
drought treatment and subsequent recovery fully sug- 
gested that the seedlings from Pulan had stronger re- 
sistance to the same degree of drought stress than 
those from Cuoqin. That was to say, S. purpurea from 
more arid area had stronger resistivity to drought 
stress; therefore, the results implied that the in- 
traspecific difference of response to drought presen- 
ted genetic effects. 

Low leaf RWC after drought caused obvious 
withering and chlorsis in plants from Cuoqin that suf- 
fered from greater mortality as a result, indicating 
that plants from Cuogin had been exceeded physio- 
logical drought limit. Water shortage can induce 
stomatal closure , which will influence photosynthesis 
process directly. The photosynthetic apparatus asso- 
ciated with PS II is very sensitive to drought and va- 
rious other stresses ( Yordanov et al., 1999; Mi- 
hailova et al., 2011) , which usually reduces chloro- 
phyll fluorescence. In the present study, the chloro- 
phyll fluorescence (F/F,,) of plants from Pulan 
dropped less and recovered more than those from 
Cuoqin, indicating that plants from Cuogin were 
more susceptible to drought and more difficult to re- 
cover from the stress. The decline of F/F „ often re- 
sulted from a non-radiative process of thermally dis- 


sipating absorbed light, which played a central role 


in leaf photoprotection under drought ( Chaves et 
al., 2003). If excess energy from absorbed light was 
not coped with by either photosynthesis or photore- 
spiration and the thermal dissipation, the production 
of highly reactive molecules, generated within the 
chloroplast, could be exacerbated to cause oxidative 
damage to the pthotosynthetic apparatus ( Smirnoff, 
1998; Niyogi, 1999). 

Oxidative stress is a state of damage caused by 
ROS, which consists of H,O,, 0, , OH and 'O, 
(Chaves et al., 2003). Under normal conditions, 
ROS are scavenged by antioxidant molecules and en- 
zymes that are located in different cell compart- 
ments, which maintain ROS homeostasis in plants 
(Li et al., 2014). In this work, we tested the ROS 
(H,0, & O,”) accumulation in conjunction with the 
antioxidant enzyme activities (SOD, CAT & APX) 
and the expression of their corresponding regulatory 
genes. Smaller amounts of H,O, and O,” were ob- 
served both after drought treatment and subsequent 
recovery in plants from Pulan, the reason could be 
explained by the higher activities of antioxidant en- 
zymes. For the plants from Cuoqin , because the low- 
er antioxidant enzymes activities had poor efficiency 
to remove the excessive ROS, so the plants failed to 
maintain the normal physiological state and even 
died heavily due to the oxidative damage partly. In- 
terestingly , the changes of antioxidant enzyme activi- 
ties were very similar to the expression results of the 
corresponding genes, implying that the divergences 
between plants from two sites were internal genetic 
properties. When the ROS dynamic equilibrium is 
broken, a series of oxidation products will be gener- 
ated, such as MDA, a commonly used index that 
reflected grades of cellular oxidation (Li et al., 
2014). In the present study, less content of MDA in 
the seedlings of S. purpurea from Pulan also indirect- 
ly illustrated their better capability to cope with 
drought stress. 

Osmotic adjustment is considered a crucial process 
in plant adaptation to drought because of its role in 


sustaining tissue metabolic activity and enabling re- 
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growth upon rewetting ( Morgan, 1984). Proline, 
the synthesis of which is mainly regulated by P5CS 
gene (Hare et al., 1999) , is one of the most studied 
compatible solutes in osmotic adjustment. Proline 
helps to stabilize macromolecules, protect enzymes, 
and store carbon and nitrogen for use during stress 
regimes in plants ( Ashraf and Foolad, 2007). Our 
results showed that the proline production could be 
rapidly and substantially induced under drought 
stress in S. purpurea, and the accumulation amount 
was much more in plants from Pulan, which had 
better drought tolerance. Interestingly, the change 
trend of proline content was consistent with the ex- 
pression of P5CS gene in the process. This suggested 
that extremely efficient start of regulatory genes, in- 
cluding P5CS, and accumulation of proline had sig- 
nificant roles in response and resistance to drought in 
S. purpurea. 

Transcription factors (TFs) have been demon- 
strated to play important roles at various levels in the 
signaling web to enable plants to deal with water- 
stress/drought (Tripathi et al., 2014). WRKY TFs, 
both positive and negative regulators of gene expression 
( Eulgem and Somssich , 2007) , have been reported to 
play pivotal roles in regulating many stress reactions 
in plants (Chen et al., 2012; Rushton et al., 2012). 
Specially , overexpression of several WRKY TFs has 
been proved to enhance drought and some other 
stresses tolerance in different species (Zhou et al., 
2008; Qiu and Yu, 2009; Wu et al., 2009; Moon 
et al., 2014). In the study, three WRKY TFs were 
all induced under drought stress and had higher ex- 
pression in S. purpurea from Pulan rather than Cuo- 
qin, implying their important functions in response 
to drought stress in S. purpurea and the role of grea- 
ter expression in strengthening drought tolerance in 
plants from Pulan. DREBs, which control stress-in- 
ducible gene expression in the abscisic acid ( ABA) - 
independent pathway ( Yamaguchi-Shinozaki and Shin- 
ozaki, 2006), are another class of TFs that have 
been reported to play an important role in plant re- 


sponse to water stress in many species (Liu et al., 


2013). The expression changes of two DREB TFs in 
our study also indicated their significant functions in 
response to drought stress and the role of greater ex- 
pression in driving stronger drought tolerance in 
S. purpurea from Pulan. 

The regulation of aquaporins is one important 
response of plant cells to water stress ( Luu et al., 
2007) , which tightly control the transcellular water 
movement by their amount and activity ( Chaumont 
and Tyerman, 2014). PIPs are one of the five sub- 
groups of aquaporins (Johanson et al., 2001; Daniel- 
son and Johanson, 2008) , which can be divided in- 
to two major groups in plants, PIP1 and PIP2, 
based on their sequences and water-channel activity 
(Kelly et al., 2014). Sufficient evidence has dem- 
onstrated that overexpression of PIPs significantly 
enhanced the rates of growth, transpiration and pto- 
tosynthesis in Arabidopsis and tomato plants ( Aha- 
ron et al., 2003; Flexas et al., 2006; Sade et al., 
2010), and the antisense suppression usually had 
the opposite effects (Siefriz et al., 2002; Uehlein et 
al., 2003; Siefritz et al., 2004). In addition, a pu- 
tative aquaporin gene V{/PIP1, isolated from Vicia 
faba leaf epidermis, was found to be induced in ex- 
pression by ABA and polyethylene glycol 6 000 and 
its expression might improve drought resistance of 
the transgenic plants (Cui et al., 2008). Like the 
previous studies, two PIPs genes, PIP1 and PIP2, 
were both induced in expression in S. purpurea under 
drought stress; furthermore, their expression were 
much higher in plants from Pulan, which displayed 
greater drought tolerance in the study. The results 
once again indicated that PIPs could improve the 
drought resistance of plants. 

Above all, our study indicated that different S. pu- 
rpurea populations had variant responses to drought 
stress throughout long-term adaptation and evolution. 
Rapid and efficient start of response and protection 
mechanism at physiological and molecular levels 
might be the foundation of stronger resistance to 
drought stress for S. purpurea. The results indicated 


the differential responses to climate change among 
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different S. purpurea populations. The study helps us 
to understand more about the adaptation and evolu- 
tion of alpine plants in the natural habitat on the Ti- 


betan Plateau under global climate change. 
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